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AUGMENTED BURST-ERROR CORRECTION 
FOR UNICON LASER MEMORY 

Raymond S. Lim 

Institute for Advanced Computation 
Ames Research Center, NASA, Moffett Field, Calif. 94035 

SUMMARY 

This paper describes a proposed augmented single-burst-error cor- 
rection system for data stored in the UNICON Laser Memory. In the 
proposed system, a long Fire Code with code length n > 16,768 bits is 
used as an outer code to augment an existing inner shorter Fire Code 
for burst-error corrections. The inner Fire Code is a (80,64) code 
shortened from the (630,614) code, and it is used to correct a single- 
burst error on a per- word basis with burst length b 5 6. The outer 
code, with b < 12, would be usert to correct a single-burst error on a 
per-page basis, where a page consists of 512 32-bit words. In the pro- 
posed system, the encoding and error detection processes would be 
implemented by hardware. A minicomputer, currently used as a UNICON 
memory management processor, would be used on a time-demanding basis for 
error correction. Based upon existing error statistics, this combina- 
tion of an inner code and an outer code would enable the UNICON system 
to obtain a very low error rate in spite of flaws affecting the recorded 
data. The approach of the long Fire Code described here is also 
applicable to other mass memory systems (such as rotating disks) where 
single-burst-error correction can be obtained at very low redundancy. 



I 

i -i . - 1 


INTRODUCTION 

The UNICON laser memory, or UNICON, is a digital memory manufactured 
by Precision Instrument Co. for the ILI.IAC-IV computer as part of the 
mass I/O storage. The UNICON js an on-line storage capacity of about 
10* ‘ bits. In such a high-density data storage system, contamination 
and other defects can easily obliterate a group of daca bits. In order 
to operate successfully in spit of this problem in the UNICON, an 
elaborate error-correction syst m is used. 

Error statistics obtained on the operation and checkout of the 
UNICON indicate that the desired net error rate o*' one bad word of data 
per 2 x 10** bits read might met if two Fire Codes were used for 
single-burst-error correction. These two Fire Codes would be arranged 
as a system consisting of an inner code and an outer code. The inner 
code would be a short (80,64) Fire Code used to correct an error burst 
of length b i 6 on a per-word basis. The outer code would be a long 
(16803,16768) Fire Code used to correct an error burst of length b i 12 
on a per-page basis. It should be noted that this system of an inner 
code and an outer code would not function in the same sense as the 
concepts of concatenated codes (ref. 1). Furthermore, the Rced- 
Solomon code from C«F(2 m ) successfully used in the IBM Photo-Digital 
memory (ref. 2) is not feasible for the UNICON because of the long 
codeword length of n=16803. 

The inner (80,64) Fire Code currently in use is implemented all in 
hardware by shift registers and error correction is done on-line as 
data are read off the strip. The (16803, l'>768) cuter Fire Code would 
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he decoded in a manner suggested by Chien (ref. 3). A similar scheme 
for bunt error correction with b < 11 was implemented successfully in 
the IBM model 3830 controller for the Model 3330 disk. The problem of 
effectively implementing the (16803,16768) outer code has been solved 
by a number o*' innovations and in a hybrid manner. Most important is 
the use of hardware for encoding, while using a minicomputer on a time- 
demanding basis for error correction. The minicomputer performs a 
calculation to find the burst error location i by applying two well- 
known theorems from Number Theory: the Euclidean algorithm for 

division and the Chinese remainder theorem for simultaneous congruences, 
sometimes thought to be me. . 1 y of academic interest only to Number 
Theorists. 


CODING REQUIREMENT FOR UN I CON 

The basic storage element of the UN1C0N is a metal coated plastic 
strip about 4.75 x 31 in. The recording density across the strip is 
11805 tracks over 3.8 in., or about 3000 tracks per in. Each track 
can store up to eight records of data, with each record containing 
16,768 data bits plus other record identification bits. The method for 
recording data on the strip is tu use the laser to either burn a hole 
for a 1-bit, or not to burn a hole for a 0-bit. 

Data entering and leaving the UNICON are handled by a UNICON 
Controller (ref. 4). The Controller has two interfaces: a central 

memory (CM) interface for data, and a minicomputer interface for 
control. The minicomputer (a PDP-11) functions as the UNICON memory 
management processor (UMP). The UMP obtains its control programs 
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through its own virtua' hardware interface into the CM. By means of 
this victual address space, UMP can address up to 409b pages in CM. 

Data transfer to the UNICON is always one page of data at a time from 
CM, where a page of data consists of 512WX32B, rfhich is 16,384 hits. 

For each page of data, the controller attaches 16WX16B of header in- 
formation and 8WX16B of checksum parity bits to form a data record of 
16,768 bits. The checksum is simply a single-parity check with check 
symbols from GF(2*'**). 

The basic word length in the UNICON is 64 bits. As previously 
mentiored, the UNICON currently has a built-in (80,64) Fire Code which 
is used to correct a single-burst error with burst length b < 6 in a 
per-word-basis. Because of optic imperfections and other mechanical 
flaws, data read back from the strip always contain errors. When the 
UNICON is properly aligned, most errors are correctable by the (80,64) 
Fire Code. Preliminary observations on the error statistic indicate 
that the probability of getting an error not correctable t s Fire 
Code is about one for every 10** bits of data transferred, wi.ich is about 
one error per strip. The observed uncorrectable errors, for those that 
are understood and the cause can be explained, are typically single- 
burst errors with b i 11. Thus if a long Fire Code with b i 12 is used 
to augment the existing (80,64) Fire Code to correct a single-burst- 
error on a per-page basis, then it is possible that the error rate can 
be reduced to a very low level. The long Fire Code chosen is a 
(16803,16758) code. In this arrangement, the (16803,16768) code is 
called the outer code, while the (80,64) code is called the inner code. 


r f 


DESCRIPTION OF THE FIRE CODE 

An error burst of length b is defined as a vector whose non-zero 
components are confined to b consecutive bit positions, the first and 
last of which are non-zero. It is known that cyclic codes for single- 
burst-error correction can be systematically constructed (refs. S and 6). 
These codes are rather efficient and they can be implemented with 
feedback shift registers. An important class of these codes is known as 
the Fire Codes (ref. 7). In the discussion which follows, let 
n » length of codeword in bits, 

k ■ number of message bits in a codeword, 

n-k ■ number of redundant bits in a codeword. 

b ■ maximum length of single error burst which can be corrected, 

in bits. 

m(x) • message word. 

v(x) • transmitted codeword. 

E(x) ■ error burst. 

R(x) ■ v(x) ♦ E(x) « received codeword. 

It is clear that for a given k and b, the objective is to construct 
an (n,k) code with as small a redundancy n-k as possible. 

It is also known from the Reiger bound (ref. 8) that the number of 
parity check bits of a b-burst-error correcting code is 

n-k 1 i'b (1) 

This is an upper bound, and codes which meet the Reiger bound are said 
to be optima 1 . The ratio 
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is a measure of the burst-correcting efficiency of the code. An 
optimal code has Z>1. The fire code at best has 


2b „ 2 


1 ' 

which is not very optimal with respect to the Keiger bound. 


(3) 


The Fire Cod » is a linear cyclic co-'e which can he systematically 
constructed for correcting a single Uirst of error in a codeword of n 
bits. The constructed codeword is in systematic form; that is. the 
n-k parity check bits are simply concatenated after the message bits. 

A Fire Code with code symbols from GF (2) which is capable of 
coirecting any burst of length b or less and of simultaneously 
detecting any burst of length d - b is best described by its generator 
polynomial 


g(x) - p (x) lx ♦ 1) 


(4) 


where 


p(x) * an irreducible polynomial of degree m whose roots have 
order e; that is, the period of p(x) is e. 

c 1 b ♦ d ■ 1 

d Zb 

m a b 

(c.e) * 1, that is, c and e are relatively prime. 

For pure error correction purpose, the best choice is to set m=b and 
d=b. This results in a Fire Code with the following parameters: 
n ■ LCM (e,c) « ec 
n-k *c+m*3b-l 
k « ec - (3b - 1) 
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Another version of the Fire Code, which is the ha’f- length version, has 
the following code paruments: 

£ 

n ■ (j) e, where c is an even integer, 
n-k - 3b - 2 

k • (f) e - (3b - 2) 

* 

In the conventional decoder implemented by shift registers, 
decoding the Fire Code requires two n shifts, n shifts for parity 
checking and another n shifts for error correction. As the first 
n shifts attri^-ited to parity checking are concurrent with the reading 
r ration, no time delay occurs in this Deration. The second n 
shifts needed to locate the error burst and correct it is a time delay 
due to error correction. If n is very large, like 16768 bits, then 
this long time delay may be intolerable in real life operation. 

However, if there exists within easy reach some power for general 
computation, resided either within the host computer or within the 
controller of the storage system, then the second n shifts can be 
reduced to the minimum of (e«-c-2) plus computation time. 

THE (80,64) FIRE CODE AND ITS IMPLEMENTATION 
The generator polynomial for the (80,64) Fire Code As 

g(x) « (x 6 ♦ x ♦ 1) (x 10 ♦ 1) (5) 

In this cis*, p(x) - (x ft ♦ x ^1) is primitive so the order of 

its roots is 
* * 

e * 2 b - 1 * 63 (6) 
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The code parameters for this code are as follows: 

(63,10) • 1, that is, 63 and 10 are relatively prime, 
n • ec ■ 63 x 10 • 630 
n-k ■c*m»10«’6»l6 
k ■ 630 - 16 - 614 
m ■ 6 so that b ■ 6 
c ■ 10 ■ b ♦ d • 1 

■ 6 ♦ d - 1 so that d • 5 

Note that the condition d £ b is not met. This means that this code 
probably will not be able to correct all possible single-burst errors 
of length b i 6. 

From the above calculations, the code generated by g(x) in 
equation (5) is a (630,614) code. Since the word length of the UNIC0N 
is 64 bits, the code must be shortened to a (8C,64) code by making the 
(614 - 64) • 550 most significant bits of the message equal to zero and 
omitting them. This process will not affect the encoding and parity 
checking calculation since leading zeros will not affect them. It does, 
however, affect the error correction procedure since the unaltered 
procedure would require 550 shifts corresponding to the 550 omitted 
zeros before reading the actual received codeword out of the buffer. 
Instead of making the buffer wait for 550 shifts, the standard 
procedure is to premultiply the decoding shift register by 
modulo g(x). 

The parity-check calculation, as it stands, is the residue of 
X lb R(x) modulo g(x). An additional automatic multiplication by 
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x 550 is desired; that is, the residue X 5bb is desired. With the 

help of a program written in basic, the remainder r(x) after 

dividing X 5bb by g(x) is found to be 

, . , 13 12 11 6 3 2 , 

r(x) »(x *x 4 x ♦x^xtx^x) 


This basic program is shown in figure 1. The encoding shift register 
is shown in figure 2, and the decoding shift register is shown in 
figure 3. 


THE (16803, 16761 > r IRE CODE AND ITS IMPLEMENTATION 

The objective of this code is to have k ■ 16768 and b ■ 12. 

The procedure to find a Fire Code that meets this requirement is as 
fol lows: 

1. Let m ■ d ■ b ■ 12 and c ■ 2b - 1 » 23 
so that c^m«23+12»35 

2. g(x) - p(x) (x 23 ♦ 1) 

, n k ♦ (n-k) 16803 C1 , A 

c c i 3 

Thus p(x) has degree m ■ 12, and a period e i 514 is 

12 5 

Jesired. With some looking, the trinomial (X ♦ X 4l) is 
not primitive but has a period e « 819 is acceptable. Thus 
g(x) - (X 12 ♦ X 5 4 1) (X 23 4 1). 

Now that g(x) is found, the code parameters are: 
n ■ e* * 819 x 23 ■ 18837 
n-k = c ♦ m ■ 23 ♦ li ■ 35 


k * (18837 - 35) « 18802 


The result is a (18837,18802) code, and this code can he shortened to 
(16803,16768) by (witting 2034 leading zeros. 

The error correction part of the decoding requires n • 16803 
shifts, which is too long a tine delay for practical application. 
However, there exists a high-speed decoder for this code suggested by 
Chien (ref. 3). This high-speed decoder consists of two feedback shift 
registers, one based on the factor (X ♦ 1 ) and the other based on the 
factor p(x). The two registers are run in synchronism for parity 
checking and error detection. A general diagram of this decoder is 
shown in figure 4. In this decoder, the syndrome is presented by the 
remainder of (X c »l) and p(x), and the error pattcr.i is identified 
when it has the same representation in both registers. The operation 
of this decoder is as follows: 

1. Shift both registers n times to enter the received 
codeword R(x). 

2. If R(x) has no error, the syndrome in both registers is 
zero. 

3. If the syndrome is not zero, shift only the (X^ ♦ 1 ) register 
until its (b- 1 ) high-order bit positions equal to zero. 

The error pattern is contained in the b low-order bit 
positions. The maximum number of shifts is (c-1). Freeze 
this register. 

4. Next shift the p(x) register until its contain is matched 
with the error pattern in the (X c ♦ 1 ) register. The maximum 
number of shifts is (e-1). If a correctable error occurs, 
this match condition will occur, otherwise a" uncorrectable 


error is detected. 


I 


’ T T 1 r * 
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From the above description, it is noted that the maximum decoding time 
is equal to the minimum of (•♦c-2), as compared to n for the con- 
ventional decoder. 

Suppose a correctable burst of length * b at location 


i, o < i < n, has occurred, then the error can be written as 

E(x) ■ x 1 B(x) 

The remainders calculated in the registers are 


S p (x) - x 1 B(x) 
S c (x) - x‘ B(x) 


modulo p(x) 
modulo (x" ♦ 1 ) 


(*) 

(9a) 

(9b) 


It is no^ed that the syndrome S^(x) is capable of determining the 
burst pattern and its location up to a multiple of c. Therefore, by 
shifting the syndrome S^(x) in its register and testing for zero at 
the (b-1) high-order bit positions, the error pattern and its location 
i can be determined quickly up to multiple of c. That is, 


i = r 


modulo c 


(10a) 


where o i r^ < c. Now that the error pattern is determined, the p(x) 

register can be shifted until a match is obtained. This gives 

i 5 r mouulo e (10b) 

P 

If i can be solved in equation (10), then the decoding is complete. 
The problem in equation (10) is a classical problem of simultaneous 
congruence in Number Theory. If e and c are relatively prime, then 
the Chinese Remainder Theorem can be used to solve for i. 

The Chinese Remainder Theorem simply states that 
"Numbers which satisfy the simultaneous congruence i - a ] 

(mod vij),. . ., i = a^ (mod p^)exist if the are relatively 


- 11 - 


prime in pairs, and such numbers constitute a single number 

j 

class modulo * W k ." 

1 

Now to solve for i. According to the Euclidean division algorithm. 


if (e,c) ■ 1, there exists intergers 

A and A such that 
C p 


A e ♦ A c ■ 1 
c p 

modulo (n ■ ec) 

(ID 

Multiply both sides by i, the result 

is 


1 ■ (A. •) i ♦ (A p C) i 

modulo n 

(12) 

From the shift registers, note that 

Q 

i ■ q c ♦ r for (x ♦ 

x c 

1) register 

(13a) 

i ■ q p e ♦ r p for p(x) 

register 

(13b) 

By substituting equation (13) into equation (12), the result is 


i * (A e) r ♦ (A c) r 
v c / c P P 

modulo n 

(14) 

Since 

A q e c = 0 modulo n 

c X 


(15a) 

A q_ e c r 0 modulo n 

P P 


(15b) 


The key equation is equation (14). The computation of and A p can 

be done easily off line by the UMP with the numbers A e and A c 
' ' c p 

stored in its memory. Once the error-pattern-match condition is found 
in the decoder, and r p are available and the UMP can compute i 

hy equation (14). Once i is solved, UMP will fetch the appropriate 
words from CM and perform the error correction with the error pattern 
stored in the p(x) register. 

It is not convenient to describe the detailed implementation of the 
(16803,16768) Fire Code since its code length is too long, 16803 bits. 
However, in order to illustrate the principle of this high-speed decoding 
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scheme, a shorter (35,27) Fire Coda with b S 3 shortened to a i 24 ,16) 
code is fully illustrated in the Appendix. 

CONCLUSION 

This paper has presented an augmented scheme for single-burst-error 
correction for the UNICON laser memory. The augmentation is the 
addition of an outer code, a (16803,16768) Fire Code, to correct for a 
single error burst with length b 1 12 on a per-page basis. This 
(16803,16768) Fire Code is an addition to the existing (80,64) Fire 
Code which has already been implemented in the UNICON to correct a 
single error burst with length b 5 6 in a per-word-basis. Based upon 
observed error statistics, this augmented scheme would reduce the 
error rate to about one err.r per 2 x 10 11 bits of data transfer, which 
is about >ne error per two strips read. The (16803,16768) code would 
be implemented by a hybrid method using a high-speed hardware decoder 
for syndrome calculation and a minicomputer for error location 
computation and error correction. For the case where the periods of 
the ti a polynomial factors of g(x) are relatively prime, a simple 
application of the Chinese Remainder Theorem would determine the burst 
locat ion. 
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APPENDIX 


The purpose of this appendix is to illustrate the principle of the 
high-speed decoder by implementing a (24,16) Fire Code shortened from 
the (35,27) code. The generator polynomial for this code is: 
g(x) - (x 3 ♦ x ♦ 1) (x 5 ♦ 1) 

where p(x) « (x 3 ♦ x ■ 1) is primitive and hence its order e ■ 7. 

The parameters of the code are: 

in ■ b ■ 3 , d ■ b ■ 3 


c ■ 2b - 1 * 5 
n « ec “ 7 x 5 » 35 
n-k«c + m»3b-l»8 
k - n - (n-k) ■ 27 


The computation of A^e 

and A c 
P 

are as follows: 

A e ♦ A 

c ■ 1 

mod 35 

c p 



7 A ♦ 5A 

> 1 

mod 35 


c p 


A =3 and A * 3 
c p 

Thus 


i * (21 r £ ♦ 15 Tp) mod 35 

Lets shorten the (35,27) code to a (24,16) code by making the leading 
11 message bits equal to zero. Encoding is not affected. However, 
decoding must account for the missing 11 leading zeros. Any multiple 
of C * 5 shifts of the (x ♦ 1) register has no effect. Any 
multiple of e ■ 7 shifts of the p (x) register has no effect. 

For the (x c ♦ 1) register, 


(3 ♦ 11) - (5 x 2) • 4 shifts 
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meaning input data should be connected to the inprt of the 5th stage 
of the (x c ♦ 1) register. For the p(x) register, 

(3 ♦ 11) - (7 x 2) « 0 shifts 

meaning input data should bt connected to the input of the first stage 
of the shift register. The decoder is shown in figure A-l. 

For a message m(x) of 



0 

123 

456 

789 10 11 12 

13 

14 15 

m(x) ■ 

1 

011 

000 

101 101 

0 

0 1 

the encoded 

V (x) 

is 






0 

123 

456 

789 

10 

11 

12 

13 

14 

15 

1 

1 

01 

234 

567 

L 

011 

000 

101 

1 

' 1 "1“ 

0 

1 

0 

0 

1 

i: 

1" 

Oil 
1 

101 


m(x) n-k 


In these examples, bit-0 is the MSB and it is sent first. To illustrate 
the principle of error correction, lets introduce two error patterns as 
follows: 

Error Pattern No. 1 

0 123 9 10 11 12 13 15 I 01 7 

I 

1 1 1 0 .... 0 | 00. .. .0 


Ej(x) 


0 000 0 


10 


Error Pattern No. 2 
0 123 

E 2 (x) « 0 000 


4 

0 


56 7 


11 1 


89 

00 


15 I 01 7 

I 

0 1 00. ...0 

I 


i ■ 5 

The corrections for these two error patterns are shown in figures 
A-2 and A-3, respectively. 
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01001 pem--ihis program compute? x(expn) moo g(x> by linfar shift 

01002 HfcM-.Bj CISTFB SEQUENCES, 

2101P R£M-. INITIALIZATION 

01029 PRINT "N"j 

0X0J0 INPUT N 

01040 LET |»1 

01045 DIM *(60) 

01050 r UP J»0 TO 16 

01091 X(J)«0 

01052 r.C XT J 

01060 MEM--9HIFT RIGHT 9UBT, 

01065 X(16)«X(1") 

01066 *ab)>*(14) 

01067 *(14)»*(13) 

01068 X(13)RX(12) 

0106V X(i?)*X(11) 

01060 IF *(lfr)«X(i0) GOTO 63 

01061 >(ll)'l 

01062 GOTO 90 

01083 *(m«c 

01090 IF *(]6)sX(9) GOTO 93 

01091 X(10)*l 

01092 GOTO 100 

01093 X ( 1 0 ) *0 

01100 X(9)«X(S) 

CJllO X(R)*X(7) 

"1120 X(7)aX(f>) 

01130 IF X(lfc)«X(S) GOTO 133 

01131 X ( 6 ) ■ 1 

01132 GOTO 140 

01133 XCftJRl* 

01140 X(5)«X(4) 

01190 X(4)»*(3) 

01160 X(3)RX(2) 

Gil 70 *(?)«X(1) 

01180 IF X(1M«X(0) GOTO 183 

01181 X(l)*l 

01182 GOTO 190 

01190 IF x ( 1 6 ) ■ I goto 193 ^RBCEDINO PAGE BLANK NOT FILMED 

01191 X (0) * 1 

C1192 GOTO 21« 

01193 X(U)i0 

01200 REM— END OF SHIFT RIGHT SUBT. 

01205 LFT X ( 1 6 ) * 0 
01210 LET t»0 
01220 LET N»N«i 
01230 IF N>«0 GOTO 6b 

01240 PRINT '’XCH), X<1), X ( 1 5 ) " 

01241 PRINT X(0)|X(l)fX(2)|X(3)|X(4)|X(5)|X(6)|X(7) 

01242 PRINT X(8)»X(O}|X(10)|X(ll)»X(l?)|X(11)|X(14)»X(15) 

01250 PRINT 

01260 GOTO 20 

P1999 END 


Figure 1.- A basic program to calculate X modulo 
g(x)“X 16 *X* 1 *X*°+’X 6 *X+l . 
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and premultiplication r(r)*(X 



Figure 4.- High-speed decoder for the fire code. 








decoder for the (24,16) fire code. 




















Figure A - 2 . 2 .- p(x) decoding for error pattern J 
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Figure A-5.2.- p(x) decoding for error pattern 






